BABAR-CONF-02/012 
SLAC-PUB-9282 
July 2002 



Measurement of the Inclusive Electron Spectrum in Charmless 
Semileptonic B Decays Near the Kinematic Endpoint 



The BABAR Collaboration 
July 24, 2002 



Abstract 

A preliminary result on the study of the inclusive electron spectrum in S — >• Xucv decays above 
the kinematic limit for the dominant B XcCv transitions is presented. This study is performed 
at the PEP-II asymmetric B Factory, where B meson pairs are produced in the decay of the T{4S) 
resonance. For the electron momentum range of 2.3 — 2.6 GeV/c in the T{4S) rest frame, the 
partial branching ratio is measured to be AB{B Xueu) = (0.152 it 0.014 it 0.014) • 10"'^. 
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1 Introduction 



The principal physics goal of the BABAR experiment is to establish CP-violation in B mesons and 
to test whether the observed effects are consistent with the prediction of the Standard Model. In 
this model CP violating effects are predicted from the CKM matrix of the couplings of the charged 
weak current to quarks. A precise determination of the matrix element \Vub\ will place constraints 
on the unitarity of the CKM matrix and thus the consistency with the minimal Standard Model. 

The extraction of \ Vub\ is a challenge, both theoretically and experimentally. While at the parton 
level weak interactions can be reliably calculated, meson decays depend on the b quark mass and 
its motion inside the B meson. Theoretical calculations of the semileptonic decay rate in terms 
of \Vub\ rely on an operator product expansion (OPE) in inverse powers of the b quark mass ||l| 
and thus depend on the choice of renormalization scale and include non-perturbative parameters. 
Experimentally, the principal difficulty is the separation of -B ^ Xucv decays from the dominant 
B —I- XcCu decays. Selection criteria applied to achieve this separation generally make it difficult 
to translate the observed rate to the full decay rate. 

In this paper we present a measurement of the inclusive electron spectrum for charmless semilep- 
tonic B decays in the momentum range of 2.3 — 2.6 GeV/c as measured in the T(4S') rest frame. In 
the rest frame of the B meson, the kinematic endpoint of the electron spectrum for the dominant 
B X^ev decays is ~ 2.3GeV/c and ~ 2.6GeV/c for B X^ev decays. In the rest frame of 
the T{4:S), the B mesons have a momentum of ~ 0.3GeV/c and thus the electron spectrum is 
convoluted by a spread of ib0.2GeV/c, extending the endpoint to higher energies. Nevertheless, 
a narrow interval of about 300MeV/c remains that is dominated by electrons from B — > Xueu 
transitions. This interval covers approximately 10% of the total electron spectrum for charmless 
semileptonic B decays. The extrapolation from the limited momentum range near the endpoint to 
the full spectrum is a very difficult task because the OPE breaks down in this part of phase space. 

This analysis is based on the same method as previous measurements of the lepton spectrum 
near the endpoint |, | (ARGUS), |, |] (CLEO). 

2 Detector and Data Sample 

This analysis is based on data recorded in 1999-2000 with the BABAR detector at the PEP-II energy 
asymmetric e"^e^ collider at the Stanford Linear Accelerator Center. The data sample corresponds 
to an integrated luminosity of 20.6 fb~^ that was collected at the T(4S) resonance (ON), plus an 
additional sample of 2.6 fb~^ that was recorded about 40 MeV below the T(4S) peak (OFF). 

The relative normalization of the ON and OFF-peak data is 7.87 it 0.01 it 0.04. This factor 
has been derived from luminosity measurements, which are based on the QED cross section for 
e"'"e~ fJ''^fJ'~ production, corrected for the energy dependence of the cross section. The systematic 
error on the relative normalization is estimated to be ~ 0.5%. This error accounts for small 
variations in the detector response and is significantly smaller than the systematic error of ~ 1.6% 
on the absolute fJ-'^fJ-^ cross section measurement. 

The BABAR detector has been described in detail elsewhere Q. The most important components 
for this study are the charged particle tracking system, consisting of a five-layer silicon detector 
and a 40-layer drift chamber, and the electromagnetic calorimeter assembled from 6580 CsI(Tl) 
crystals. Electron candidates are selected on the basis of the ratio of the energy detected in the 
calorimeter to the track momentum, the calorimeter shower shape, the energy loss in the drift 
chamber, and the angle reconstructed in the ring imaging Cherenkov detector. 
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The electron identification efficiency and the probabilities to misidentify a pion, kaon, or proton 
as an electron have been measured with clean samples of tracks that were selected from the data. 
This experimental information is introduced into the Monte Carlo simulation to improve the agree- 
ment with the data. Tracking efficiencies and resolution have been studied. A comparison with 
the simulation has revealed small differences, which have been taken into account. No significant 
impact of non-Gaussian resolution tails has been found in the endpoint region. 

The Monte Carlo simulation of B ^ X^ev events is based on the ISGW2 model Q. In the 
current version, the hadrons are represented by single particles or resonances with masses up 
to 1.5 GeV/c^ and non-resonant contributions are not included. For B — > Xcev transitions three 
models are employed to simulate different decay modes. The decay to D*eu is modeled following a 
form factor based parameterization of HQET for decays to Dev and higher mass charm meson 
states the ISWG2 model is used. The non-resonant decays to D^*^TTeu are modeled according to a 
prescription by Goity and Roberts @. 



3 Data Analysis 
3.1 Event Selection 

For this analysis, electron candidates are selected in the momentum range from 1.5 to 3.5 GeV/c 
in the T{4S) rest frame. The solid angle is restricted by the electromagnetic calorimeter coverage, 
defined by the laboratory polar angle range of —0.72 < cos^iab < 0.92. In this momentum and 
angular range, the efficiency for identifying an electron has been measured to be epiD = 0.91 ±0.02, 
while the average hadron misidentification probability is less than 0.2%. The selected electron 
sample is dominated by electrons from semileptonic decays of B and D mesons, non-resonant qq 
production and QED processes. In addition, photon conversions and Dalitz decays contribute 
background at low momenta and J/tp decays contribute at higher momenta. Furthermore, there 
are sizable contributions from hadrons misidentified as electrons. To suppress low- multiplicity 
QED processes, including r+r~ pairs, the number of charged tracks per event is required to be 
greater than three. This background and non-resonant hadronic events are further suppressed by 



a restriction on the ratio of Fox-Wolfram moments H2/HQ < 0.4 |10|. In semileptonic decays, the 
neutrino carries a substantial momentum. In events in which the only undetected particle is this 
neutrino, its direction can be inferred from the missing momentum in the event, defined as the 
difference between the net momentum of the two colliding beam particles and the vector sum of 
all detected particles, charged and neutral. Therefore, the selection of these decays can be greatly 
enhanced by requiring that the measured missing momentum exceed 1 GeV/c and point into the 
detector fiducial volume. Furthermore, since the B mesons are produced almost at rest and the 
high momentum electron and neutrino point in nearly opposite directions, we require that the 
angle between the electron candidate and the missing momentum be greater than 7r/2. Candidate 
electrons are rejected if, when paired with an opposite-sign electron, the invariant mass of the pair 
is consistent with J/ip mass, 3.05 < Mg+g- < 3.15 GeV/c^. 

The detection efficiencies are estimated using Monte Carlo simulation. For the selection cri- 
teria described above, the detection efficiency for charmless semileptonic decays in the electron 
momentum interval of 1.5 — 2.7 GeV/c ranges from ~ 0.4 to ~ 0.25. 
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3.2 Background Subtraction 

The raw spectrum of the highest momentum electron in events selected by the criteria described 
above is shown in Figure |l]a, separately for data recorded ON and OFF resonance. For the OFF- 
resonance data, the momenta are scaled to compensate for the 0.4% difference in the c.m.s. energies 
of the two data samples. 

Also shown in Figure ^ is a fit to approximate the non-resonant background contribution. For 
the purpose of this fit the normalized OFF-peak data in the interval of Pcms = 1-5 — 3.5 GeV/c are 
combined with the ON-peak data in the interval of Pcms = 2.7 — 3.5 GeV/c above the kinematic 
limit for B XuGv decays. The x^/dof for a fit with a 4-th degree Chebyshev polynomial is 
51/51. Figure |l]b shows the residuals of the fitted function from the unsealed OFF-peak data. The 
fit significantly reduces the statistical uncertainty in the estimate of the continuum background 
that is to be subtracted. 

The result of the subtraction of the fitted continuum background is shown in Figure Hi. Also 
shown are the Monte Carlo predictions of the expected signal from B X^ev decays and back- 
ground contributions from all other processes, B X„ei^. In this simulation, the B — > Xueif 
branching ratio is assumed to be 1 • 10^^ and the relative normalization of data and Monte Carlo 
simulation is set by the number of electrons in the momentum interval between 1.5 and 2.3 GeV/c. 
The result of the subtraction of all backgrounds is shown in Figure ^3. 

The number of selected events, split into three momentum intervals, together with the estimated 
background contributions are listed in Table |l[ The intervals are chosen to emphasize differences 
in the background composition: the interval 2.0 and 2.3 GeV/c is dominated by semileptonic B 
XcCu decays, with a small contribution from hadron misidentification, the 2.3 — 2.6 GeV/c interval 
is mostly populated by -B — > X^ev decays and continuum background, while the 2.7 — 3.0 GeV/c 
interval contains almost exclusively continuum background. 

Table 1: The number of signal and background events and signal selection efficiencies e{B X^ev) 
for three different intervals in the electron momentum. A^oN refers to the number of selected 
electrons recorded on the T(4S') resonance, A^qff is the fitted number of continuum background 
events. The N{B — > XcCv) and N{B — > Xc — > e) refer to the Monte Carlo estimate of the remain- 
ing semileptonic background, from prompt and secondary decays (including ^/)(2S')), respectively. 
N{B —5- J /tp e^e~) is an estimate of J/ip background remaining after the mass cut. N{B Xc) 
is a background from hadronic B decays to charm. N{B X^ev) is the resulting number of signal 
electrons. 



M omentum pe ( GeV/c) 


2.0-2.3 


2.3-2.6 


2.7-3.0 




74, 140 ± 272 


6, 455 ±80 


1,932 ±44 


iVoFF 


7, 749 ± 165 


4, 051 ±93 


1,903 ±37 


NiB Xcev) 


61, 158 ±470 


470 ± 41 





N{B J/i> ^ e+e") 


666 ± 49 


128 ±22 





N{B Xc^e) 


338 ± 35 


18 ±8 





N{B Xc) - mis-ID 


373 ± 37 


92 ±18 


4±4 


N{B Xucv) 


3, 857 ±572 


1, 696 ± 133 


25 ±57 


e{B^Xuey){%) 


33.9 ±1.0 


27.7 ± 1.3 


28.6 ±22.9 
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Figure 1: Electron momentum spectrum in the T{4S) rest frame: (a) ON-peak data (solid circles), 
OFF-peak data scaled by the ON/OFF-ratio (open circles), (b) OFF-peak data (unsealed) after 
the subtraction of the fitted continuum. The line show the result of a fit to the continuum spectrum 
(using both ON- and OFF-peak data) in the interval Pcms = 1-5 — 3.5 GeV/c with a 4-th degree 
Chebyshev polynomial. 



3.3 Determination of S — > XuCv Branching Ratio 

For a given interval in the electron momentum, the inclusive partial branching ratio is calculated 
according to 

^^^ iV0N-iV0FF-iVB^X.e. ^,^^^^)_ (1) 

Here Nq^^ refers to the number of electrons detected ON-peak and A^off refers to the fitted contin- 
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Figure 2: The eleetron momentum spectrum in the T{4S) rest frame: (a) ON-peak data after con- 
tinuum subtraction (soUd triangles), and MC predicted background from BB events {B X^eu) 
(open triangles), (b) ON peak data after subtraction of continuum and Monte Carlo predicted 
(B Xy^ev) backgrounds (data points with statistical errors). For comparison, the histograms 
show the expected signal spectrum from B — > Xucv decays. 



uum background in a specified momentum interval, Nj^^Xueu 

is the background from BB events 
derived from Monte Carlo simulation normalized to the total number of electrons in the momentum 
range 1.5 — 2.3 GeV/c, e is the total efficiency for detecting a signal electron from B X^eu decays 
(including bremsstrahlung in detector material) , and S^ad accounts for the distortion of the electron 
spectrum due to final-state radiation. This last is a momentum dependent correction, derived from 
the Monte Carlo simulation, amounting to ~ 11% in the range 2.3 — 2.6 GeV/c. As the overall 
normalization the total number of produced BB events is used, N^^ = (22, 630 it 19 ± 362) • 10^. 
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The differential branching ratio as a function of the electron momentum in the T(45') rest frame is 
shown in Figure ^. Integrated over the interval from 2.3 to 2.6 GeV/c, the partial branching ratio 
is (statistical error only): 

AB{B Xucv) = (0.152 ± 0.014) • 10~^ (2) 



0.1 



^ 0.075 
> 

o 



m 
I 

O 



Oh 



0.05 



0.025 







BaBar 



2.1 2.2 2.3 2.4 2.5 2.6 2.7 

Electron Momentum (GeV/c) 



Figure 3: The branching ratio B{B X^eu) as a function of the electron momentum in the T(45') 
rest frame. The data (statistical errors only) are compared to a prediction (solid line) based on the 
ISGW2 model assuming a total inclusive branching ratio of 1 • 10~^ for B Xueu decays with X^ 
masses up to 1.5 GeV/c^. The spectrum is corrected for final-state radiation and bremsstrahlung. 
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3.4 Systematic Errors 

The following sources of systematic errors have been considered: 

• the Monte Carlo simulation of the signal events and the impact on the detection efficiency, 

• continuum background subtraction, 

• the subtraction of all other backgrounds based on Monte Carlo simulation of the production 
process and the detector simulation, 

• the uncertainty of the B meson momentum spectrum in the T(4S') rest frame. 

The systematic error introduced by the efficiency estimation for the signal events is caused 
by uncertainties in simulation of event selection (4%), the track reconstruction (1%), and the 
electron identification efficiencies (2%). The error on the selection efficiency depends not only on 
the accuracy of the detector model, but also on the adequacy of the ISGW2 model on which the 
event generation is based. The latter error has been estimated from the impact of the variation 
(within their current uncertainties) of individual branching ratios for B meson decays and variation 
of the event selection cuts. The total systematic error on the efficiency amounts to 5%. 

An error in the continuum background subtraction can be introduced by the choice of the fitting 
function and the fitting procedure. The impact of the choice of the function has been estimated 
by trying different functional forms, such as an exponential function and Chebyshev polynomials. 
Possible differences in continuum spectra in ON-peak and OFF-peak data have been assessed by 
including or excluding ON-peak data above B X^ev kinematic limit in the fit. The observed 
variations in the number of continuum events leads us to attribute a 2% systematic uncertainty on 
the number of continuum events. This is equivalent to 5% systematic error in the number of signal 
events. 

The 6 — > c background consists of electrons and misidentified hadrons (see Table ||) . The domi- 
nant contribution are electrons from semileptonic B XcCu decays. There is also a much smaller 
contribution of secondary electrons from B Xc Xgev and B — s- ip{2S) ee decays, but very 
few with momenta above 2.3 GeV/c. The error in the estimate of prompt lepton background arises 
primarily from the uncertainty in the decay dynamics and the branching fractions of the individual 
decay modes. The impact of these uncertainties (4%) is estimated by variation of these branching 
ratios within their current uncertainties. In addition, the uncertainty in the electron identification 
(2%) has been included. The systematic error (28%) on the background from misidentified hadrons 
arises from both the uncertainty in the spectrum of charged hadrons above 2.3 GeV/c and the 
overall uncertainty on the mis-identification probabilities for pions and kaons. In addition, there 
is a small correction for electrons from J/-^ decays that are not vetoed by the mass requirement; 
the uncertainty in this correction is 20%. Background from photon conversions and Dalitz decays 
have been analyzed and found to be negligible above 1.5 GeV/c. The total systematic error on the 
Monte Carlo based 6 — > c background subtraction is estimated to be 6%, translating to a relative 
error of 3% on the number of signal events. 

Since the T{4S) mass is so close to the threshold for BB production, the B meson momentum 
and thereby the endpoint of the electron spectrum from B X^ev decays is highly sensitive to 
the total energy of the colliding beams. Thus variations of the colliding beam energy introduces a 
systematic error in the B — > Xcev background subtraction. These variations are included in the 
Monte Carlo simulations, and the impact of the difference between data and Monte Carlo is studied 
using a sample of fully reconstructed B mesons. The comparison of the measured and simulated 
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momentum distributions, currently limited by statistics, shows that the systematic error from this 
effect does not exceed 5%. 

In summary, Table ^ contains list of systematic errors for the momentum range from 2.3 to 
2.6 GeV/c; the total systematic error on the partial branching ratio measurement for this momen- 
tum range is ~ 9%. 

Table 2: Systematic error on the AB for the momentum range from 2.3 to 2.6 GeV/c. All numbers 
are given in percent. 



contribution 


Sx/x 


S{AB)/AB 


efBciency e{B XuCiy) 


5 


5 


continuum subtraction 


2 


5 


background B Xcev 


4 


1 


background B J/tp — > e+e^ 


20 


2 


background B Xc - mis-ID 


28 


2 


B movement 


5 


5 


Nbb ^ 


1.6 


1.6 


radiative corrections 


10 


1 


total 
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4 Conclusion 

The result of this analysis, the fully corrected differential branching ratio as a function of the 
electron momentum in the T(4S') rest frame, is presented in Figure |3[ Integrating over the interval 
from 2.3 to 2.6 GeV/c results in the partial branching ratio 

ABiB Xuev) = (0.152 ± 0.014 ± 0.014) • 10"^ (3) 

This result agrees very well with the measurement by the CLEO collaboration The measured 
electron spectrum in the endpoint region is well reproduced by the ISGW2 model. 

5 Extraction of ll^^l based on CLEO Measurements 

To determine the charmless semileptonic branching fraction B{B X^ev) from the partial branch- 
ing fraction AB{Ap), one needs to know the fraction /„(Ap) of the spectrum that falls into the 
momentum interval Ap. Heavy Quark theory describes the Fermi-motion of the quarks inside the 
meson in terms of a shape function that depends on non-perturbative QCD. To leading order, the 
same shape function describes all h — s- qlv transitions (here q represents any light quark). 

The CLEO collaboration [|l^] has recently used the measurement of the inclusive photon spec- 
trum from 6 — > S7 transitions to derive the parameters describing the shape function and to 
calculate the lepton momentum spectrum for B — > X^ev transition. They quote a value of 
fu{Ap) = 0.074 ± 0.014 ± 0.009 for the interval Ap from 2.3 to 2.6 GeV/c. 

Relying on the CLEO measurement, the result presented here translates into a total branching 
ratio of 

B{B ^ Xuev) = (2.05 ± Q.27exp ± 0.46/J • 10-^ (4) 
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Based on studies developed independently by two groups |12]-[17|, we adopt the following rela- 
tionship for the extraction of \Vuh\ from the total branching fraction (revised most recently at 
the 2002 CKM workshop) 



Kb I = 0.00445 



B{B XJv) 1.55ps 
0.002 ~ 




(1.0 ±0.020 ±0.052) 



(5) 



where the first error arises from the uncertainty in the OPE expansion, and the second from the 
uncertainty in the h quark mass, assuming mb(l GeV) = 4.58 ± 0.09 GeV. Based on this expression 
and our measurement of the branching fraction, we find 



Here the first error is the measurement uncertainty from the combined statistical and systematic 
error, and the second refers to the uncertainty on the extraction of \Vub\ from the branching ratio 
as stated above. The remaining errors are taken from the CLEO analysis: the third refers to the 
determination of fu and the fourth represents an estimate of the validity of the assumption that the 
shape function can be extracted from the b ^ s'j spectrum. The uncertainty due to the theoretical 
assumption of quark-hadron duality remains unquantifiable. 

The BABAR collaboration is planning to develop this and other methods for extracting the total 
branching ratio and \Vub\ in the near future. 
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